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Abstract. Interstellar and atmospheric extinctions redden the observational pho- 
tometric data and they should be handled rigorously. This paper simulates the effect 
of reddening for the modest case of two main sequence Ti = 6500K and T2 = 5500K 
components of a detached eclipsing binary system. It is shown that simply subtract- 
ing a constant from measured magnitudes (the approach often used in the field of 
eclipsing binaries) to account for reddening should be avoided. Simplified treatment 
of the reddening introduces systematics that reaches ~ O.Olmag for the simulated 
case, but can be as high as ~ 0.2mag for e.g. B8 V-K4 III systems. With rigorous 
treatment, it is possible to uniquely determine the color excess value i5(B — V) from 
multi-color photometric light curves of eclipsing binaries. 
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1. Introduction 

Although interstellar extinction has been discussed in many papers 
and quantitatively determined by dedicated missions (lUE, 2MASS 
and others), there is a lack of proper handling in the field of eclipsing 
binaries. The usually adopted approach is to calculate the amount 
of reddening from the observed object's coordinates and its inferred 
distance and to subtract it uniformly, regardless of phase, from pho- 
tometric observations. This paper shows why this approach may be 
inadequate, especially for objects where interstellar extinction and the 
color difference between both components are significant. Atmospheric 
extinction is a better-posed problem: similarly as interstellar extinction 
depends on E(Ji — V), atmospheric extinction depends on air-mass, 
which is a measurable quantity, whereas £'(B — V) has to be estimated. 



2. Simulation 

To estimate the effect of reddening on eclipsing binaries, we built a syn- 
thetic binary star model, consisting of two main sequence G9 V-F5 V 
stars with Ti = 5500i^, Ri = O.861i?0 and T2 = 6500K, R2 = l.SSGi?© 
and 1 day orbital period. The simulation logic is as follows: for the given 
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Figure 1. The reddening law adopted from Cardelli et al. (1989). 

phase, we calculate the effective spectrum of the binary by convolving 
Doppler-shifted individual spectra of the visible surfaces of both com- 
ponents. To this intrinsic spectrum we rigorously apply interstellar and 
atmospheric extinctions (both as functions of wavelength). We then 
convolve this reddened spectrum with instrumental response function 
(composed of the filter transmittivity and detector response functions) 
and integrate over the bandpass wavelength range to obtain the flux. 
In contrast, we use the same intrinsic spectrum without rigorously 
applying the reddening. To simulate the subtraction^ of a reddening 
constant from photometric observations, we simply divide the intrinsic 
spectrum by the flux that corresponds to this constant. Finally, we 
calculate the flux in the same manner as before and compare it to the 
flux obtained by applying rigorous reddening. 

For building synthetic light curves we use PHOEBE^ (Prsa &; Zwitter, 
2004; in preparation). Each light curve consists of 300 points uniformly 
distributed over the whole orbital phase range. To be able to evalu- 
ate the impact of reddening on photometric light curves exclusively, 
all second-order effects (limb darkening, gravity brightening, reflection 
effect) have been turned off. 

We take Kurucz's synthetic spectra {R = 20000) from precalculated 
tables by Munari et al. (2004; in preparation). The used {UBV)j{RI)c 
response data (filter x detector) are taken from ADPS (Moro & Mu- 



^ To emphasize that this constant is subtracted from photometric observations 
(given in magnitudes), we shall use the term constant subtraction throughout this 
study, although we are dividing when operating with fluxes. 

PHOEBE stands for PHysics Of Eclipsing BinariEs and it is based on the Wilson- 
Devinney code. See http://www.f iz.uni-lj .si/phoebe for details. 
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Figure 2. Left: The change of the effective wavelength of the Jolmson B filter due 
to reddening of the simulated G9 V-F5 V binary. Right: The (B — V)j color index 
on the G9 V-F5 V (5500 K-6500 K) temperature interval, calculated by integrating 
the spectrum over both filter bandpasses. 

nari, 2000), where we apply a cubic spline fit to obtain the instrumental 
response function. 

For interstellar extinction, we use the Cardelli ct al. (1989) empirical 
formula (Fig. 1), where TZy = 3.1 was assumed throughout this study. 
Schlegel et al. (1998) interstellar dust catalog was used to obtain the 
maximum color excess £^(B — V) values for different lines of sight. 

For atmospheric extinction we use the equation triplet for Rayleigh- 
ozone-aerosol extrinction sources given by Forbes et al. (1996) and 
summarized by Pakstiene Sz Solheim (2003). The observatory altitude 
h = km and the zenith air-mass are assumed throughout the study. 

To rigorously deredden the observations for the given TZy and E(B — 
V), it is necessary to determine the reddening for each wavelength of 
the spectrum. Correcting differentially and integrating over the filter 
bandpass then yields the dereddcned flux of the given filter. However, 
without spectral observations, it is difficult to calculate properly the 
flux correction. Since Cardelli et al. (1989) formula depends on the 
wavelength, the usually adopted approach found in literature is to 
use the effective wavelength Agg of the filter transmittivity curve to 
calculate the reddening correction. We demonstrate the implications in 
the following section. 

3. Implications 

3.1. Interstellar extinction 

By comparing the rigorously calculated fluxes against intrinsic fluxes 
with a simple constant subtracted, we come to the following conclu- 
sions: 1) taking the effective wavelength of the filter bandpass should 
be avoided. Since the flux is the integral over the filter bandpass, X^s 
has a conceptually different meaning. Furthermore, Aeff of the given 
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Figure 3. Left: The discrepancy due to simplified constant subtraction approach 
(solid line) compared to the rigorously applied reddening (points) for a G9 V-F5 V 
binary. The subtracted constant was obtained from the effective wavelength 
(-^off ~ 4410. SA) of the Johnson B transmittivity curve. Right: Overplotted light 
curves with the subtraction constant calculated so that the magnitudes in quarter 
phase are aligned. There is still a measurable difference in eclipse depth of both light 
curves. i5(B — V) = 1 is assumed. 



filter depends heavily on the effective temperature of the observed 
object and on the color excess i?(B — V) (Fig. 2). To determine the 
subtraction constant, one has to make sure that the integral (rather 
than any particular wavelength) of the both curves is the same. Fig. 3 
shows the discrepancy between the properly calculated light curve and 
the one obtained by subtracting a Agfr-calculated constant. Table I 
summarizes the differences between the proper treatment and other 
approaches. 2) Even if the subtraction constant is properly calculated, 
the light curves still exhibit measurable differences in both minima 
(Figs. 3, 4). This is due to the effective temperature change of the 
binary system during eclipses. For the analysed case, the difference in 
B magnitude is ~ O.Olmag, which is generally observable. 3) If light 
curves in three or more photometric filters are available, it is possible 
to uniquely determine the color excess value -E'(B — V) by comparing 
different color indices in-and-out of eclipse. The reddening may thus 




Figure 4- The difference between the rigorously calculated reddening and a constant 
subtraction in Johnson B filter for the G9 V-F5 V binary (left) and a B8 V-K4 III 
binary (right) during primary minimum. i?(B — V) = 1 is assumed. 
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Tabic I. The summary of different approaches to calculate the wave- 
length to be used for the dereddening constant. Xes is determined by 
requiring that the area under the spectrum on both sides is equal. 
Ams is the value of extinction in B filter and EAm^ is the devia- 
tion from the rigorously calculated value. All values are calculated 
for E(B — V) = 1 at quarter phase. Note that Ams is smaller than 
TZv+E(B—Y) = 4.1, since our simulated binary is cooler than lOOOOK. 



Approach: 


Acfi 






Rigorously calculated value: 


4470.2A 


4.03 


0.00 


Filter transmittivity: 


4410.8A 


4.09 


0.06 


Filter transmittivity -|- reddening law: 


4452.1 A 


4.05 


0.02 


Intrinsic spectrum: 


4436. 3 A 


4.06 


0.03 


Effective (reddened) spectrum: 


4583.6A 


3.90 


-0.13 



be properly introduced to the fitting scheme of the eclipsing binary- 
analysis program. This was done in PHOEBE. 

3.2. Atmospheric extinction 

Atmospheric extinction is comprised of three different sources: the 
Rayleigh scattering, the aerosol scattering and ozone absorption (Forbes 
et al., 1996). It depends on the wavelength of the observed light and 
on the air-mass of obsevations, which is of course the same for both 
binary components; we are thus left only with wavelength dependence 
at some inferred air-mass. To assess the impact on the photometric 
data, we compare the flux change imposed on the intrinsic spectrum 
by reddening and by atmospheric extinction (Fig. 5) . We conclude that 
for weakly and moderately reddened eclipsing binaries the atmospheric 
extinction dominates the blue parts of the spectrum due to Rayleigh 
scattering, but for larger color excesses (£^(B — V) > 0.5) the reddening 
is dominant throughout the spectrum. 



4. Discussion 

By not properly taking reddening into account, we introduce system- 

atics of ~ O.Olmag into the solution for our simulated binary. One 
may ask: is this difference worth bothering with? To demonstrate that 
it is, let us simulate a bit more exotic eclipsing binary with blue B8 
main-sequence (Ti = 12000K, Ri = 3.221i20) and red K4 type III 
giant (T2 = 4000K, R2 = 30. Oi?©) components with orbital period 
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Figure 5. The percentage of the flux change imposed on the intrinsic spectrum by 
reddening (thick line) and atmospheric extinction (thin line). If the color excess is 
large, the reddening completely dominates the spectrum (left), but if the color excess 
is moderate, the blue part of the spectrum is dominated by the Rayleigh scattering 
(right). 

of 15 days. Fig. 4 (right) shows that the discrepancy between the 
rigorously calculated light curve and the constant-subtracted one is 
as large as ~ 0.2mag, which means ~ 10% error in the determined 
distance to the observed object. On the other hand, the reddening effect 
may be reduced by using narrower filter-sets, e.g. Stromgren ubvy set. 
Furthermore, typical observed color excesses rarely exceed few tenths 
of the magnitude, which additionally diminishes this effect. However, 
future space scanning missions such as GAIA (Ferryman et al., 2001) will 
acquire measurements of eclipsing binaries without spatial bias. Thus 
interstellar extinction should be carefully and rigorously introduced 
into the reduction pipeline for eclipsing binaries. 
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